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Abstract 

The ANTARES telescope has the capability to detect neutrinos produced in as- 
trophysical transient sources. Potential sources include gamma-ray bursts, core col- 
lapse supernovae, and flaring active galactic nuclei. To enhance the sensitivity of 
ANTARES to such sources, a new detection method based on coincident obser- 
vations of neutrinos and optical signals has been developed. A fast online muon 
track reconstruction is used to trigger a network of small automatic optical tele- 
scopes. Such alerts are generated for special events, such as two or more neutrinos, 
coincident in time and direction, or single neutrinos of very high energy. 

Key words: ANTARES, Neutrino astronomy, Transient sources, Optical follow-up 
PACS 95.55.Vj 



1 Introduction 



The detection of high energy cosmic neutrino from a source would be direct 
evidence of the presence of hadronic acceleration within the source and provide 
important information on the origin of the high energy cosmic rays. Powerful 
sources of transient nature, such as gamma ray bursts or core collapse super- 
novae, offer one of the most promising perspectives for the detection of cosmic 
neutrinos as, due to their short duration, they are essentially background free. 
For example, several authors predict the emission of neutrinos in correlation 
with mult i- wavelength signals, e.g. the Fireball model of GRBs [1]. As neu- 
trino telescopes observe a full hemisphere of the sky (even the whole sky if 
downgoing events are considered) at all times, they are particularly well suited 
for the detection of transient phenomena. 

In this paper, the implementation of a strategy for the detection of transient 
sources is presented. This method, earlier proposed in [2], is based on the 

^ Also at University of Leiden, the Netherlands 
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optical follow-up of selected neutrino events very shortly after their detection 
(Section 3) by the ANTARES neutrino telescope [3]. The alert system, known 
as "TAToO" (Telescopes and ANTARES Target of Opportunity) [4], uses an 
online track reconstruction with a pointing accuracy of about 0.5 degrees. This 
reconstruction algorithm is described in Section 4. Its characteristics allow the 
triggering of optical telescopes such as TAROT [5] and ROTSE [6] . In order to 
improve the precision of these alerts, an additional reconstruction algorithm, 
described in Section 5, which takes into account the detailed detector geometry 
is used offline. 



2 The ANTARES detector and data acquisition 

The ANTARES neutrino telescope is located in the Mediterranean Sea, 40 km 
from the coast of Toulon, France, at a depth of 2475 m. The detector is an 
array of photomultiplier tubes (PMTs) arranged on 12 slender detection lines, 
anchored to the sea bed and kept taut by a buoy at the top. Each line comprises 
up to 25 storeys of triplets of optical modules (OMs), each housing a single 
10" PMT. Since lines are subject to the sea current and can change shape 
and orientation, a positioning system comprising hydrophones and compass- 
tiltmeters is used to monitor the detector geometry. Data taking started in 
2006 with the operation of the first line of the detector. The construction of 
the 12 line detector was completed in May 2008. The main goal of the experi- 
ment is to search for neutrinos of astrophysical origin by detecting high energy 
muons (>100 GeV) induced by their neutrino charged current interaction in 
the vicinity of the detector. Due to the large background from downgoing cos- 
mic ray induced muons, the detector is optimised for the detection of upgoing 
neutrino induced muon tracks. 

2.1 Data acquisition 

The task of the ANTARES data acquisition system (DAQ) [7] is to collect 
the data from all the individual PMTs of the detector and pass them to the 
filtering algorithms which search for a collection of signals compatible with a 
muon track crossing the detector. 

From the DAQ point of view, each storey is an independent acquisition unit in- 
cluding a processor, buffering RAM and an Ethernet link to the shore station. 
Individual PMT pulses above a threshold of typically 0.3 photoelectron (re- 
ferred to as 'hits') are digitized offshore in the form of 'hit time' and 'charge', 
and sent to a computer farm onshore for further processing. Due to biolumi- 
nescence activity and 40 K decays in the sea water, each PMT has an average 
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counting rate of the order of 100 kHz, requiring large bandwidth for data trans- 
mission to shore. The data are transmitted through gigabit Ethernet links on 
single optical fibers to the computing farm where the filtering algorithms are 
executed, reducing the event rate to a few tens of hertz. The overall system 
is supervised by a state machine which handles the various commands needed 
to configure, start and stop the acquisition on both the offshore and onshore 
processors. 

Data are time structured in the form of time slices of 104.85 ms, allowing the 
data of the full detector for the duration of one time slice to be sent to a 
single computing node. The synchronization of the 300 offshore processors is 
performed by a 20 MHz clock distribution system broadcasted to all storeys. 
In particular, the start of any data taking period is stamped using an external 
GPS signal giving the absolute timing at the location of the detector, allowing 
an absolute time accuracy better than 1 //s [8]. 

2.2 Data filtering 

The main goal of the filtering algorithm is to select hits compatible with 
the propagation of Cherenkov light emitted by a muon crossing the detector, 
among the background from bioluminescence and 40 K decays characterized 
by uncorrelated hits on single OMs. The filtering is based on local 'clusters', 
defined either as coincidence hits on OMs of the same storey within a narrow 
time window or as a single hit with a large amplitude. There are two main 
filtering algorithms running simultaneously, both searching for a combination 
of local clusters within a typical 2.2 fis time window. The first algorithm 
requires five local causally connected clusters anywhere in the detector, while 
the second requires at least two local clusters in adjacent or next-to-adjacent 
storeys. All hits within a few microseconds around these clusters define an 
"event" and are kept for further online and offline reconstructions. 

2.3 The TAToO trigger 

Once a single data filtering computer has processed its time slice, the result- 
ing events are sent to a data distribution service to which the "Data Writer" 
storing task, monitoring tasks and the online event reconstruction [9] are con- 
nected. The "alert" application described in Section 4 analyses the data stream 
of reconstructed events, selects candidates fulfilling various criteria, and then 
generates the TAToO alerts. 

An important performance parameter for the alert system is the time between 
the crossing of the detector by a high energy muon and the time at which 
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an alert is sent. This time is the sum of the data dispatching time from the 
offshore photomultipliers to the onshore computing farm (1.5 s) and of the data 
processing time of an entire time slice by the filtering algorithm (<5 s). For the 
time-being, an extra buffering time of 60 seconds is added; this delay will be 
removed after an upgrade of the acquisition to improve the performance of the 
alert sending capability. The time needed to reconstruct the event direction 
and verify the alert criteria typically amounts to a few milliseconds. Therefore, 
the total delay between an interesting particle crossing the detector and a 
TAToO alert is currently 66 s. 



3 The TAToO Run Control 

The TAToO Run Control (RC) is a stand-alone Qt control application!"^! which 
channels the triggers generated by the alert application to the optical telescope 
network. 

The connections to this network are checked periodically and automatic re- 
connection is performed resulting in a fully autonomous and stable system. A 
veto prevents an alert to be sent if the ANTARES event counting rate exceeds 
a given threshold. In addition, if the alert criteria are fulfilled soon after a 
previous alert has already been issued, the new alert is stored in a FIFO and 
sent only after a certain period of time. This time lag, currently set at one 
hour, is used to avoid alert pileup in the optical telescope network. Manual 
alerts can also be generated and sent. All alerts are sent using the Gamma-ray 
bursts Coordinates Network (GCN) [10] normalized format, allowing easy im- 
plementation of connections to additional telescopes. Information about the 
event that triggered the alert, i.e. a unique identifier, the time and the celes- 
tial coordinates, the number of hits used in the reconstruction and the track 
reconstruction quality are sent to the optical telescopes network at the time 
of the alert. 



4 The alert criteria 

The criteria for the TAToO trigger are based on the features of the neutrino 
signal from the expected sources. Several models predict the production of 
high energy neutrinos greater than 1 TeV from GRBs [11] and from Core 
Collapse Supernovae [12]. Under certain conditions, multiplet of neutrinos 
can be expected [13]. 



5 Software framework originally from Trolltech, now Nokia. 
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Two online neutrino trigger criteria are currently implemented in the TAToO 
alert system: 

• the detection of at least two neutrino induced muons coming from similar 
directions within a predefined time window; 

• the detection of a single high energy neutrino induced muon. 

A basic requirement for the coincident observation of a neutrino and an optical 
counterpart is that the pointing accuracy of the neutrino telescope should be 
at least comparable to the field of view of the optical telescope. 



4-1 The online track reconstruction algorithm 

To select the events which might trigger an alert, a fast and robust algorithm 
is used online to reconstruct tracks from the calibrated data. This algorithm 
uses an idealized detector geometry which does not rely on the dynamical 
positioning alignment. As a result, the hits of the three OMs of a storey 
are grouped and their location assigned to the barycenter of the storey. The 
storey orientations as well as the line-shape deviations from straight lines 
are not considered in the online reconstruction. A detailed description of this 
algorithm and its performance is found in Ref. [9]. The principle is to minimize 
a x 2 which compares the times of selected hits with the expectation from a 
Cherenkov signal generated by a muon track. The resulting direction of the 
reconstructed muon track is available within 10 ms and the obtained minimum 
X 2 is used as a fit quality parameter to remove badly reconstructed tracks. 



4-2 Neutrino selection criteria 

Atmospheric muons, whose abundance at the ANTARES detector [14] is roughly 
six orders of magnitude larger than the one of muons induced by atmospheric 
neutrinos, are the main background for the alerts and have to be efficiently sup- 
pressed. Among the surviving events, neutrino candidates with an increased 
probability to be of cosmic origin are selected [15]. 

4-2.1 Atmospheric muon background rejection 

Atmospheric muons resulting from the interaction of cosmic rays with nuclei 
in the atmosphere represent the main component of the background. Atmo- 
spheric muons propagate downgoing through the detector and can be sup- 
pressed with an elevation cut, selecting only the upgoing events. However, some 
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badly reconstructed atmospheric muons classified as upgoing may remain, and 
quality cuts are applied to reduce this contamination to an acceptable level. 

In order to establish the criteria used for our neutrino selection, we have anal- 
ysed a subsample of data taken by ANTARES after the completion of the 
12-line detector, corresponding to a livetime of 70.3 days. During this period, 
around 350 upgoing neutrino candidates were reconstructed and have been 
compared to a Monte Carlo (MC) simulation of atmospheric muons and neu- 
trinos using the same livetime. Downgoing atmospheric muons were simulated 
with Corsika [16], and normalized to match the data. The primary particle 
flux was composed of several nuclei according to Ref. [17] and the QGSJET 
hadronic model [18] was used for the shower development. Upgoing neutrinos 
were simulated according to the parameterization of the atmospheric neutrino 
flux from Ref. [19]. Only charged current interactions of neutrinos and an- 
tineutrinos were considered. The Cherenkov light produced in the vicinity of 
the ANTARES detector was propagated taking into account light absorption 
and scattering in sea water [20]. The angular acceptance, quantum efficiency 
and other characteristics of the optical modules were taken from Ref. [21]. 




Fig. 1. Track fit quality (x 2 ) distribution for all upgoing events reconstructed on at 
least 2 lines. The atmospheric muon Monte Carlo distribution has been rescaled to 
match the data. 

Figure 1 shows the distribution of the track fit quality (x 2 ) of the minimiza- 
tion procedure for all upgoing events reconstructed with at least two lines. A 
cut on the track fit quality is applied to reduce the number of atmospheric 
muons reconstructed as upgoing in the final sample. Because the fit quality 
is correlated to the number of hits used in the fit, the selection cut on the fit 
quality parameter is set to a different value according to the number of hits 
used to reconstruct the event: \ 2 < 1.3 + [0.04 * (N hit — 5)] 2 . 

Figure 2 shows the elevation distributions both for data and simulated at- 
mospheric neutrino and muon samples, after the cut on the track fit quality. 
Atmospheric muons reconstructed as upgoing are efficiently rejected and a 
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neutrino purity better than 90% is achieved. 




Fig. 2. Elevation angle distribution after a selection cut on the fit quality for all 
events reconstructed on at least 2 lines. 



Figure 3 shows the angular resolution of the online algorithm as a function of 
the neutrino energy for events reconstructed with different number of lines. 
This resolution is defined as the median of the space angular difference be- 
tween the incoming neutrino and the reconstructed neutrino-induced muon. 
For neutrinos with an energy higher than a few tens of TeV, an angular res- 
olution of 0.4 degree is achieved, despite of the approximations related to the 
detector geometry. For example, the inclination of the ANTARES line for a 
typical sea current of 5 cm/s induces a systematic angular deviation of less 
than 0.2 degree. 




Fig. 3. Angular resolution as a function of neutrino energy for events with tracks 
reconstructed with different number of lines. 
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4-2.2 Multi-neutrino trigger 



The typical signature of the transient emission of high energy neutrinos is a 
neutrino burst, i.e. a multiplet of neutrino events originating from the source 
in a short time window. A trigger for this event type is implemented as the 
detection of two upgoing events reconstructed with at least two lines in a 
15 minutes time window with a maximum angular difference of 3°. The time 
window was optimized to include most predictions of the neutrino emission 
by various models for transient sources. The 3° angular window was selected 
to match the convolution of the track reconstruction angular resolution and 
the field of view of the robotic optical telescopes (~ 2° x 2°). The accidental 
coincidence rate due to background events, from two uncorrelated upgoing 
atmospheric neutrinos, is estimated to be 7 x 10~ 3 coincidences per year with 
the full ANTARES detector. With such a small background, the detection of a 
doublet (triplet) in ANTARES would have a significance of about 3 (5) sigma. 

4-2.3 High energy event trigger 

Since the neutrino energy spectrum for signal events is expected to be harder 
than for atmospheric neutrinos, a cut on the reconstructed energy efficiently 
reduces the atmospheric neutrino background while most of the signal events 
are kept. 




Number of storeys 

Fig. 4. Distribution of number of storeys having at least one hit used in the recon- 
struction, for data and Monte Carlo upgoing events reconstructed on at least three 
lines. The vertical line indicates the alert selection criterium. 

The selection of the alert candidates is based on two simple energy estimators: 
the number of storeys used in the track fit and the total amplitude (in pho- 
toelectrons) of the hits in the storeys. Figures 4 and 5 show the distributions 
of the number of storeys and the amplitude, respectively, both for data and 
Monte Carlo samples. 
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The event selection for the high energy trigger has been tuned on atmospheric 
neutrinos in order to obtain a false alarm rate of about 25 alerts per year. This 
rate was agreed between ANTARES and the optical telescope collaborations. 
A requirement of at least 20 storeys on at least three lines and an amplitude 
greater than 180 photoelectrons will select 25.7 ± 0.4 high energy events per 
year with the full 12 line configuration of the ANTARES detector. 




Fig. 5. Total amplitude distribution in photoelectrons (pe) for data and Monte Carlo 
upgoing events reconstructed on at least three lines. The vertical line indicates the 
alert selection criterium. 



4.2.4 Trigger performance 

The performance of these two triggers has been studied using a neutrino Monte 
Carlo generated with an E~ 2 energy spectrum. The TAToO alert criteria select 
neutrinos of energies above 1 TeV for the multi-neutrino trigger, and above 
10 TeV for the single high energy trigger. 

The pointing accuracy of the neutrino telescope has two contributions: the 
angle between the incoming neutrino and the resulting muon, and the angle 
between the muon and the reconstructed trajectory. The former contribution 
is due to kinematics: the more energetic the neutrino, the more co-linear the 
resulting muon will tend to be with the neutrino direction. The latter con- 
tribution is determined by the performance of the track reconstruction algo- 
rithm. The bidimensional point spread function expected for neutrino events 
selected by the single high energy trigger is shown in Figure 6. As an illustra- 
tion, ~80% of the events are reconstructed within the field of view (FOV) of 
a typical robotic telescope. Figure 7 shows the dependence of this fraction on 
the neutrino energy. 
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Fig. 6. Bidimensional point spread function. The black solid-line square corresponds 
to a typical robotic telescope field of view (« 2° x 2°). 




log i0 [E v (GeV)] 



Fig. 7. Fraction of events inside the field of view of a typical robotic telescope (2° x 2°) 
as a function of the energy of the event, assuming that the tracks originated from 
the center of the field of view, for both online and offline reconstruction algorithms. 

5 Determination of the refined direction 



Although the pointing accuracy of the online reconstruction algorithm is suit- 
able for the field of view of the telescopes used for the follow-up, the use of 
the detailed knowledge of the detector geometry can further improve the de- 
termination of the neutrino direction. For this purpose, we use the standard 
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ANTARES offline reconstruction algorithm [22]. 

Since the ANTARES lines are not rigid structures, sea currents can move the 
top buoy by several metres, and distort the line positions from a vertical line 
geometry, thus affecting the direction of the reconstructed muon trajectory. 

In order to achieve the best track reconstruction performance, it is necessary 
to monitor the relative positions of all OMs with an accuracy of better than 
20 cm, equivalent to 1 ns timing precision. In addition, a precise absolute 
orientation of the whole detector is necessary to point to individual neutrino 
sources in the sky. Two independent monitoring systems are used to attain 
the required accuracy: 

• A high frequency long baseline acoustic system giving the 3D position of 
hydrophones placed along the line. These positions are obtained by trian- 
gulation from emitters anchored at the bases of the lines. 

• A set of tiltmeter-compass sensors giving the local tilt angles of each OM 
storey with respect to the horizontal plane (pitch and roll) as well as its 
orientation with respect to the Earth's magnetic field (heading). 

In order to obtain a quasi-online precise detector geometry (within a delay of 
typically few tens of minutes), the shape of the detector lines is derived from 
a model which estimates the mechanical behaviour under the influence of the 
sea water flow, obtained from online measurements of the sea current. The 
positions of the OMs are then calculated by combining the line shape with 
the measurements of the tilt and orientation angles of the storeys given by the 
tiltmeter-compass sensors. 

5. 1 The offline likelihood fit 

The offline reconstruction algorithm derives the muon track parameters that 
maximize a likelihood function built from the difference between the expected 
and the measured arrival time of the hits from the Cherenkov photons emit- 
ted along the muon track. This maximization [22] takes into account the 
Cherenkov photons that scatter in the water and the additional photons that 
are generated by secondary particles (e.g. electromagnetic showers created 
along the muon trajectory). 

The value of the log-likelihood per degree of freedom (A) from the track re- 
construction fit is a measure of the track fit quality and is used to reject badly 
reconstructed events, such as atmospheric muons that are mis-reconstructed 
as upgoing tracks. The distribution of the variable A is shown in Figure 8 for 
both data and Monte Carlo events. Only tracks reconstructed with A > —5.2 
are kept for the determination of the refined track direction. 
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Fig. 8. Track fit quality (A) distribution for upgoing events in both data and atmo- 
spheric Monte Carlo samples. 



5.2 Offline reconstruction performance 



Figure 9 shows the median Monte Carlo computed angular resolution as a 
function of the neutrino energy for events selected by the high energy alert 
trigger, for both online and offline algorithms. The improvement obtained with 
the offline algorithm is clearly visible at low energies. The angular resolution 
for neutrino energies above 10 TeV is 0.5° for the online algorithm and 0.35° 
for the offline algorithm. As illustrated in Figure 7, ~92% of the signal events 
are reconstructed within the FOV of a typical robotic telescope. 



— •— Online algorithm 




rj 

2.5 3 3.5 4 4.5 5 5.5 6 

log i0 [E v (GeV)] 



Fig. 9. Angular resolution obtained for both online and offline reconstructions as a 
function of the neutrino energy. 
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6 Optical follow-up procedure 



The optical follow-up strategy is based on short-term observations (for rapidly 
fading sources) and a long-term follow-up (mainly for core collapse SNe searches) 
Once the alert is sent at time Tq, optical images are collected as soon as possi- 
ble by the available telescopes. The fast slewing of the small robotic telescopes 
allow the collection of images 5 to 10 seconds after the reception of the alert. 
The follow-up procedure, which typically extends over 1 month with images 
taken each night close to the alert T , makes use of the refined direction. The 
observations at times later than T are organized as follows: if the refined di- 
rection is more than 0.5° away from the initial direction and has a poor offline 
reconstruction quality, the follow-up of the alert is cancelled. Otherwise, the 
pointing direction of the telescopes is updated with the refined coordinates 
and all subsequent images are centered around that direction. 



7 Summary 

The method used by the ANTARES collaboration to implement the search 
for coincidence between high energy neutrinos and transient sources followed 
by small robotic telescopes has been presented. Of particular importance for 
this alert system are the ability to reconstruct online the neutrino direction 
and to reject efficiently the background. With the described ANTARES alert 
sending capability, the connected optical telescopes can start taking images 
with a latency of the order of one minute, which will be reduced to about 15 s 
in the near future. The precision of the direction of the alert is much better 
than one degree. The quasi-online availability of a refined direction obtained 
using the measured geometry of the ANTARES detector further improves the 
quality and efficiency of the alert system. 

The alert system is operational since early 2009, and as of December 2010, 
27 alerts have been sent. After a commissioning phase in 2009, all alerts had 
an optical follow-up in 2010. These numbers are consistent with the expected 
trigger rate, after accounting for the duty cycle of the neutrino telescope. 

The optical follow-up of neutrino events significantly improves the perspective 
for the detection of transient sources. A confirmation by an optical telescope of 
a neutrino alert will not only provide information on the nature of the source 
but also improve the precision of the source direction determination in order 
to trigger other observatories (for example very large telescopes for red-shift 
measurement). The program for the follow-up of ANTARES neutrino events 
is already operational with the TAROT and ROTSE telescopes and results 
based on analysis of the optical images will be presented in a forthcoming 
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paper. This technique could be extended to observations in other wavelength 
regimes such as X-ray or radio. 
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